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Abstract. The cosmological QCD transition affects primordial density perturba- 
tions. If the QCD transition is first order, the sound speed vanishes during the 
transition and density perturbations fall freely. For scales below the Hubble radius 
at the transition the primordial Harrison-Zel'dovich spectrum of density fluctuations 
develops large peaks and dips. These peaks grow with wave number for both the 
radiation fluid and for cold dark matter (CDM). The peaks in the radiation fluid 
are wiped out during neutrino decoupling. For cold dark matter that is kinetically 
decoupled at the QCD transition (e.g. axions) these peaks lead to the formation of 
CDM clumps of masses 1O- 2O M < M < 10- 10 M Q . 



1. Introduction 

The transition from hot deconfined quarks and gluons to hot hadrons happens in the 
early Universe at a temperature T+ ~ 150 MeV [Q. Lattice QCD results for the 
physical values of the quark masses indicate that the QCD phase transition is of first 
order @. 

For a first-order QCD phase transition the sound speed vanishes during the 
coexistence of the deconfined and confined phases ||. The important scales are the 
mean distance between the hadronic bubbles at their nucleation, i? n uc, which is a 
few cm for homogeneous nucleation [Q, and the Hubble scale i?H ~ mp/I^ 2 ~ 10 4 
m. Cosmological density perturbations are affected by the QCD transition at scales 
A ~ i?H- At the coexistence temperature T* the pressure p+ = p{T±) is fixed. Therefore 
pressure gradients and the sound speed vanish for wavelengths much larger than the 
typical bubble separation. This gives rise to large peaks and dips above the primordial 
spectrum of density perturbations in the radiation fluid. These peaks grow at most 
linearly with wave number. 

At T ~ 1 MeV the neutrinos decouple from the radiation fluid. During this 
decoupling the large peaks in the radiation spectrum are wiped out by collisional 
damping §, §. 

Today the universe is dominated by dark matter, most likely cold dark matter 
(CDM). If CDM is kinetically decoupled from the radiation fluid at the QCD 
transition, the density perturbations in CDM do not suffer from the neutrino damping. 
Examples for kinetically decoupled CDM are primordial black holes, axions, or 
WIMPzillas, but not supersymmetric dark matter and heavy neutrinos 0. 
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At the time of the QCD transition the energy density of CDM is small, i.e. 
p CDM (T t ) - Kr 8 p RAD (i;). Kinetically decoupled CDM falls into the potential wells 
provided by the dominant radiation fluid. Thus, the CDM spectrum is amplified 
on subhorizon scales. The peaks in the CDM spectrum go nonlinear shortly after 
radiation-matter equality. This leads to the formation of CDM clumps with masses 
lCT 20 Af Q < M < KT 10 M Q §. 

If the QCD transition is strong enough, these clumps could be detected by 
gravitational femtolensing Q. The clumping of kinetically decoupled CDM implies 
that axions, if they are the CDM, could hide from axion search experiments ||. 

2. Kinetically decoupled cold dark matter 

Let us discuss the properties of some popular CDM candidates at the QCD scale. 

The lightest supersymmetric particle, most probably the neutralino, is one of the 
candidates for CDM jl(J. Constraints from LEP 2 and cosmology, together with the 
assumption of universality at the GUT scale, show that the neutralino mass is m x > 42 
GeV . It is essential to distinguish between the chemical freeze-out and the kinetic 
decoupling of ncutralinos. The chemical freeze-out happens when the annihilation rate 
of neutralinos drops below the Hubble rate, r ann /iJ < 1. At freeze-out the rate for 
neutralino annihilation, r ann = (va ann )n x , ^ s suppressed by the Boltzmann factor in 
the number density of the neutralinos, n x ~ (m x T) 3 / 2 exp(— m x /T). The freeze-out 
temperature of the neutralino is Tj ~ m x /20 > 2 GeV. 

Kinetic decoupling, in contrast, is determined by the elastic scattering between 
neutralinos and the dominant radiation fluid. The interaction rate for elastic scattering 
is T e i = (va e i)n, where n ~ T 3 is the number density of relativistic particles. The 
kinetic decoupling of the neutralino happens at the 10 MeV scale J5J and therefore 
neutralinos belong to the radiation fluid at the QCD transition. The same holds true 
for heavy neutrinos with < m v < 1 TcV. 

Superheavy WIMPs, WIMPzillas with masses 10 12 - 10 16 GeV @, decouple from 
the radiation fluid well before the QCD transition. Although WIMPzillas scatter with 
the radiation fluid at T*, the momentum transfer in these interactions is tiny compared 
to the WIMPzilla's momentum. Even the large number of scatterings per Hubble time 
cannot change the momenta of WIMPzillas significantly (if mwiMP > 10 6 GeV). Thus, 
WIMPzillas cannot be dragged along by the radiation fluid. 

Axions could be the dominant matter today if their mass is small, i.e. m a ~ 10~ 5 
eV, which corresponds to a breaking of the Peccei-Quinn (PQ) symmetry at the scale 
/pq ~ 10 12 GeV [Jk|. These axions could be produced coherently due to a initial 
misalignment of the axion field and by the decay of axionic strings. The initially 
misaligned axion field starts to oscillate coherently when the axion mass has grown to 
m a (Tx) ~ 3F(T X ), where T t ~ 1 GeV. Below T x ~ 1 GeV the oscillating axion field 
evolves as CDM. 

If the reheating temperature after inflation lies above /pq, the axion field is 
inhomogeneous on scales larger Rn(Ti). These inhomogeneities produce axion 
miniclusters with mass M mc ~ 1O _12 M0 and radius i? mc ~ O.li?© today Q. 

A further candidate for CDM that is kinetically decoupled at the QCD scale is 
primordial black holes (PBH) ]l5j produced before the QCD transition and therefore 
with masses A/bh *C lAf . In order to survive until today PBH should have 
Af BH > 10 15 g w ICT^AIq. PBH in the range from KT 18 M to 10~ 16 A/ Q would 
radiate too strongly to be compatible with 7-ray observations . 
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3. Sound speed 

The QCD phase transition starts with a short period (lCP 4 iH) of tiny supercooling, 
1 — Tsc/T* ~ 1CP 3 . When T reaches T sc , bubbles nucleate at mean distances i? nuc ~ 2 
cm 0| . The bubbles grow most probably by weak deflagration [|| , because the surface 
tension is very small |i"7fl . The released energy is transported into the deconfined phase 
by shock waves, which reheat the deconfined phase to within 10 -6 £h- The entropy 
production during this very short period is AS/ S ~ 1CP 6 . Further bubble nucleation 
in the reheated quark-gluon phase is prohibited after this first 10~ 4 <h- Thereafter, 
the bubbles grow adiabatically due to the expansion of the Universe (coexistence of 
the phases). The transition completes after 10 _1 tH- 

During the reversible coexistence period (isentropic) T and p are fixed, while p 
can vary, therefore 



= 0. (1) 

iscntropic 



This implies that there are no pressure gradients at scales A > i? n uc, while both phases 
coexist during the first-order phase transition. 

The thermodynamics of the cosmological QCD transition can be studied in lattice 
QCD, because the baryochemical potential is negligible in the early Universe, i.e. 
fib/T* ~ 1CP 8 . We fit the equation of state from lattice QCD |l8[. 



4. The evolution of density perturbations 

During inflation density perturbations Sp have been generated with an almost 
scale- invariant Harrison-Zel'dovich spectrum |]l9| . In the radiation dominated era 
subhorizon (fc phys ;§> H) density perturbations oscillate as acoustic waves, supported 
by the pressure Sp — <? a 8p. It is useful to study the dimensionless density contrast 
5 = Sp/ p as a function of conformal time and comoving wave number. The subhorizon 
equation of motion reads 

5" + c 2 s k 2 S = , (2) 

if we assume that the phase transition is short compared to the Hubble time JjJ . From 
the lattice QCD equation of state the period of vanishing c s lasts for O.I^h, in the bag 
model it lasts for 0.3tn- 

To illustrate the mechanism of generating large peaks above the HZ spectrum we 
discuss the subhorizon evolution of the density contrast S in the bag model. Before and 
after the QCD transition c 2 s = 1/3 and the density contrast oscillates with constant 
amplitudes A m resp. A ont . At the coexistence temperature T* the pressure gradient, 
i.e. the restoring force in Eq. (|2|), vanishes and the sound speed (0) drops to zero. 
Thus S will grow or decrease linearly, depending on the fluid velocity at the moment 
when the phase transition starts. This produces large peaks in the spectrum which 
grow linearly with the wave number, i.e. A mt /A m \ peaks = k/k%. 

Kinetically decoupled CDM falls into the gravitational potential wells of the 
radiation fluid during the coexistence regime. 

Fig. [I] shows the final spectra for the density contrast of the radiation fluid 
and of CDM for a fit Q to lattice QCD. The spectra have been obtained from 
numerical integration of the fully general relativistic equations of motion. The mass 
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Figure 1. The transfer functions for the density contrast of CDM and of the 
radiation fluid due to the QCD transition. On the horizontal axis the wave number 
k is represented by the CDM mass contained in a sphere of radius Tt/k. 

Mi corresponds to the wave number k\ of the bag model analysis. Mi coincides 
with the CDM mass inside the Hubble horizon at the QCD transition, ~ 10~ 8 M Q . 
For the equation of state from the fit to lattice QCD a WKB analysis shows that 
A out /A in \ peaks = (k/k 2 ) 3/4 , which defines the mass M 2 . We conclude that at the 
horizon scale the modification of the HZ spectrum is mild, whereas for scales much 
smaller than the horizon big amplifications are predicted. Without tilt in the COBE 
normalized spectrum the density contrast grows nonlinear for k phys /H £ 10 4 resp. 10 
for the bag model resp. lattice QCD equation of state. 

5. Clumps in CDM 

At T ~ 1 MeV the neutrinos decouple from the Hubble scale. During their decoupling 
they wash out all fluctuations of the radiation fluid below a CDM mass of 10 _5 M Q . 
Thus, at the time of nucleosynthesis all peaks in the spectrum of the radiation fluid 
(Fig. [l]) are erased and the energy density on scales A ~ -Rh(1 MeV) is homogeneous. 

For kinetically decoupled CDM collisional damping is irrelevant. Peaks in these 
types of CDM survive and grow logarithmically during the radiation era. Shortly after 
equality they grow nonlinear and collapse by gravitational virialization to clumps of 
M < 10~ 10 M©. From a COBE normalized HZ spectrum and the equation of state 
from quenched lattice QCD a clump of mass 10~ 15 M Q would have a size of ~ lOi?0. 
For a tilted spectrum with n — 1.2 the size of the same clump would be ~ IRq p). 

We find that the peaks in the CDM spectrum lead to clumps of masses 10 -20 — 
1O _1O M . Today, these clumps would have a density contrast of 10 10 — 10 17 , where 
the lower value corresponds to a 10 — 15 Mq clump from an untilted CDM spectrum, 
the bigger value is for a 1O _2O M0 clump from a tilted CDM spectrum. The evolution 
of these clumps in the late stages of structure formation remains to be investigated 
(disruption, mergers, etc.). 

For a larger amplification during the QCD transition, e.g., if it should turn out 
that the latent heat is bigger than the value from present lattice QCD calculations, 
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more compact clumps are possible. These could be subject to femto-lensing gj. With 
the values from lattice QCD, the CDM clumps are not compact enough to lie within 
the Einstein radius, which is Re ~ 0.02i? Q for a 10 _15 Mq clump. 

The clumping of CDM changes the expected rates for some dark matter searches, 
because some of the rates depend on the space-time position of the detector, star, 
or planet. Especially experiments looking for axion decay in strong magnetic fields 
would not provide a new limit on the axion mass if they find nothing. These 
experiments may just tell us that we are not sitting in an axion clump currently. 
These consequences remain to be studied further. 
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